Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 5 February 2008 



(MN M?eX style file vl.4) 



Quasi Periodic Oscillations in Low Mass X-Ray Binaries 
and Constraints on the Equation of State of Neutron Star 
Matter 



Christoph Schaab and Manfred K. Weigel 



OS 

<: 



Institut fiir theoretische Physik, Ludwig-Maximilians Universitdt Miinchen, Theresienstr. 37, D-80333 Miinchen, Germany 
email: schaab @gsm. sue. physik. uni-muenchen. de 



5 February 2008 



(N 
O 

O 



ABSTRACT 

Recently discovered quasi periodic oscillations in the X-ray brightness of low mass X- 
ray binaries are used to derive constraints on the mass of the neutron star component 
and the equation of state of neutron star matter. The observations are compared 
with models of rapidly rotating neutron stars which are calculated by means of an 
exact numerical method in full relativity. For the equations of state we select a broad 
collection of models representing different assumptions about the many-body structure 
and the complexity of the composition of super dense matter. The mass constraints 
differ from their values in the approximate treatment by ~ 10 %. Under the assumption 
that the maximum frequency of the quasi periodic oscillations originates from the 
innermost stable orbit the mass of the neutron star is in the range: M ^ 1.92 — 
2.25 Mq. Especially the quasi periodic oscillation in the Atoll-source 4U 1820-30 is 
only consistent with equations of state which are rather stiff at high densities which is 
explainable, so far, only with pure nucleonic/leptonic composition. This interpretation 
contradicts the hypothesis that the protoneutron star formed in SN 1987A collapsed to 
a black hole, since this would demand a maximum neutron star mass below 1.6 Mq. 
The recently suggested identification of quasi periodic oscillations with frequencies 
around 10 Hz with the Lense-Thirring precession of the accretion disk is found to be 
inconsistent with the models studied in this work, unless it is assumed that the first 
overtone of the precession is observed. 
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1 INTRODUCTION 

Neutron stars contain matter in one of the densest forms 
found in the universe. Their central density ranges from a 
few time the density of normal nuclear matter to about one 
order of magnitude higher, depending on the star's mass and 
the equation of state (EOS). Neutron stars therefore provide 
us with a powerful tool for exploring the properties of such 
dense matter. In the last decades, this tool was applied to, 
among other topics, the determination of the EOS of dense, 
charge neutral, /3-equilibrated matter by means of compar- 
ing the theoretical predicted properties with observations of 
neutron stars. This was attempt ed by studying, for example , 
the maximum stable star mass ( van Kerkwiik et al. , 1995 ) 



the minimum rotation period ( Friedman et al., 1986 ; Webci 
& Glendenning, 199S), or the thermal behaviour (Tsuruta 
1966 , Schaab et al. 1996 , Page 1997 ; see Balberg, Licht- 
enstadt & Cook 1998 for a recent review). 



Recently, Strohmayer et al. (1996) and Van der Klis et 
al. (1996) discovered with the Rossi X-ray Timing Explorer 
(RXTE) kilohertz quasi-periodic oscillations (QPOs) in the 
X-ray brightness of low-mass X-ray binaries (LMXBs, see 
van der Klis 1997 for a recent review). In subsequent ob- 
servations, three QPOs were often detected simultaneously 
in a given source. The frequency separation between both 
QPOs is almost constant, although the frequencies of the 
two QPOs themselves vary by several hundred Hertz. Up to 
now, the only exceptions are the Atoll sources 4U 1608-52 
and 4U 1735-44 and the Z-source Scorpius X-1 in which the 
frequen cy separation varies with the luminosity by roughly 
±15 % (iMendez et al., 1998(j ; |van der Klis fc Wijnands, 1997 



Ford et al., 1998bD. Psaltis et al. 



199S ) found that the QPO 



data of the other LMXBs are consistent both with being con- 
stant and with varying by a similar fraction as in the two 
sources quoted before. The frequency of oscillations during 
type I X-ray bursts detected in some sources are consis- 
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tent with the frequency separation of the two QPOs or with 
its first overtone. Only in the source 4U 1636-536 the aver- 
aged frequency separation, Af = 251 Hz, and the half of 
the frequency in typ e I bursts, t^burst = 581 Hz, d iffer by 
approximately 15% (Mendez & van Paradijs, 1998). In the 



power spectrum of bursts of the same source Miller ( 199J ) 
found a second signal at 290 Hz ~ l/2z/burst- These obser- 
vations strongly favour the beat frequency model where the 
frequency separation between the two QPOs originates from 
the stellar spin, whereas the higher QPO is produced by 
accreting gas in a stable, nearly circular orbit around the 
neutron star. Though, it has to be clarified how the slightly 
varying frequency separation and the small deviation of the 
frequency separation from the burst frequency can be incor- 
porated into this model. 

Beside the two kilohertz QPOs and the burst oscilla- 
tions, also QPOs with frequencies of a few tens of Hertz were 
detected in some sources. Their frequencies correlate with 
the high frequency kilohertz QPOs. These low f requency 
QPOs were interpreted by Stella & Vietri (1998) to orig- 



inate from the Lense-Thirring precession of the accreting 
disc due to th e frame dragging effect o f the rapidly spinning 
neutron star ( Lense fc Thirring, 1918 ). 

Both, the identification of the high frequency kilohertz 
QPOs with the orbital frequency of a stable circular orbit 
and the low frequency QPOs with the frame dragging fre- 
quency of the same orbit allow to constrain the mass of the 
neutro n star and also the EOS of neutron star matter (Lamb 
[Stella fc Vietri, 199^ ). If the evidence for the de- 



et al., Um 



tection of QPOs of the innermost stable circular orbit in the 
sources 4U 1608-52, 4U 1636-536 (iKaaret et al., 19971), and 



4U 1820-30 (Zhang et al., 1998b) can be confirmed by future 



observations, the constraints are rather severe, allowing only 
a few stiff EOSs. 

The frequency separation of the two kilohertz QPOs 
show that the neutron star in LMXBs are rapidly rotating 
with periods ranging from 2.5 ms to 4 ms. It can therefore be 
expected that the geometry of the neutron star and its exte- 
rior space time is non-spherical. Since the innermost stable 
orbit is located at only a few kilometres above the star's sur- 
face, the deviation from the Kerr space time are large and 
should not be neglected. In order to compare theoretical 
neutron star models with QPO observations, a completely 
general relativistic calculation of the rotating neutron star 
structure and space time geometry is therefore necessary. 

In order to study the impact and the discrimination 
power of the QPO data in greater detail, we select a broad 
collection of modern EOSs, which were obtained utilising nu- 
merous assumptions about the dynamics and composition of 
super dense matter. To mention several, these are: the many- 
body technique used to determine the EOS; the model for 
the nucleon-nucleon interaction; description of electrically 
charge neutral neutron star matter in terms of either only 
neutrons and protons in generalised chemical equilibrium (/3 
equilibrium) with electrons and muons, or nucleons, hyper- 
ons and more massive baryon states in (3 equilibrium with 
leptons; hyperon coupling strengths in matter; inclusion of 
meson (tt, K) condensation; treatment of the transition of 
confined hadronic matter into quark matter; and assump- 
tions about the true ground state of strongly interacting 
matter (i.e., absolute stability of strange quark matter rela- 
tive to baryon matter). 



The paper is organised as follows. In Sect. Q we sum- 
marise the equations which govern the space time structure 
and compare the approximate values of the orbital frequen- 
cies, the radius of the innermost stable orbit, and the Lense- 
Thirring precession frequencies with the respective values 
from the exact numerical solution of Einstein's equations. 
The physics of the EOSs is discussed in Sect. ^ The high 
frequency kilohertz QPOs and their interpretation in com- 
bination with their compatibility with the different EOSs 
are discussed in Sect. ^. The implications of the identifica- 
tion of the low frequency QPOs with Lense-Thirring preces- 
sion are presented in Sect. H. We summarise our results, the 
constraints to the neutron star masses, and the conclusions 
concerning the neutron star EOS in Sect. ^. 



2 SPACE TIME AROUND RAPIDLY 
ROTATING NEUTRON STARS 

2.1 Einstein Equations 

The stationary, axis-symmetric, and asymptotic fiat metric 
in quasi isotropic coordinates reads 



edr + e"*'(dv3 - N'^Aty + e^"(dr" + r^dr 



(1) 

where the metric coefficients g^n, = g^u{r, 9) are functions of 
r and 9 only. The metric coefficients are determined by the 
Einstein equation (c = G = 1) 

G = SttT, (2) 

and the energy-momentum conservation 

V ■ T = 0, (3) 

where T = (e-|-p)u(8)u-fpgis the stress-energy tensor of 
an ideal ffuid with the 4- velocity 



u = Fe "et + Q.Te "e^, 



(4) 



the energy density e, and the pressure p. The Lorentz factor 
r is given by u ■ u = —1, hence 

_ 1 

T ^ [l-e^'^'^-"^ {n- N'^f'^ ^ (5) 

S7 = is the angular velocity of the ffuid with respect 

to an observer at infinity. The proper velocity U o f the fluid 
with respec t to the local Eulerian Observer Do (Smarr & 
York, 1978 ) is given by the equation 



U 



gjip r 



(6) 



Note that if the fluid were at rest with respect to the Ob- 
server Do, 17 = 0, it would not necessarily be at rest for an 
inertial observer at infinity, since Q. — N"^ 7^ 0. This phe- 
nomena of dragging o f the inertial frame was first studied 
by Lense & Thirring (1918) and turns out to be important 



for the investigation of the Lense-Thirring precession of the 
accreting disc (see Sect. 5|). 

The four non-trivial Einstein equations together with 
the energy-moment um conservat ion are solved via a finite 
difference scheme ( jchaab, 1999 ) . We follow Bonazzola et 
al. (|l99^ ) in compactifying the outer space to a finite region 
by using the transformation r ^ 1/r. The boundary con- 
dition of approximate fiatness can then be exactly fulfilled. 
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The neutron star model is uniquely determined by fixing 
one of the parameters: central density, gravitational mass, 
baryon number, Kepler orbiting frequency at the star's equa- 
tor or at the marginally stable radius, as well as one of the 
parameters: angular velocity, angular momentum, or stabil- 
ity parameter t = r/|VK|. The models of maximum mass 
and/or maximum rotation velocity can also be calculated. 



2.2 Stable Circular Orbits 

Since et and are Killing vectors, the components pt and 
Pip of the 4-momentum of a freely falling particle are con- 
stants of motion and can be identified with the negative of 
the specific energy E (in units of the rest mass mo) and 
the specific angular momentum L, respectively. For a par- 
ticle motion confined to the equatorial plane the geodesic 
equation p{I = — 1 yields: 



-2u 2 
e Pr 



'V{E,L,r), 
(7) 



where V is the effec tive potential for the particle motion 
(Misner et al., 1973, pp. 655). A circular orbit, pr — 0, is 
given by the expression for the 4-velocity of the fiuid 
inside the star. Then, the specific energy E and angular 
momentum L can be expressed in terms of the Lorentz factor 
r (Eq. ^ and the proper velocity U (Eq. |^: 



The circular orbit exists if V.,- = 0, thu^ 



(8) 
(9) 

(10) 



The proper velocity [/ of a particle corotating with the star 
is given by 



2(f>^r 

The circular orbit is stable, if V^rr < with (see also Cook, 
Shapiro & Teukolksky 1992 and Datta, Thampan & Bom- 
bacci 1998 ) 



2(0-1') 



(12) 



The zero of V^rr determines the radius r™'^ of the inner- 
most stable or marginally stable circular orbit. The Kepler 
frequency of a circular orbit as measured by a distant 
observer is given by the proper velocity U through the ex- 
pression 



2-K 



"U + N'' 



and the Lense-Thirring precession frequency by 

1 



2n 



(13) 



(14) 



We use the usual convention, that -^^ represents the partial 
differential d/dx^. 



The respective values at the innermost stable orbit with ra- 
dius r™'' are denoted by v'^'^ and vff^, respectively. The cir- 
cumferential radius rdrc of the object is linked to the equa- 
torial coordinate radius r by the expression 



e r. 



(15) 



2.3 Approximate Expressions 



Though we shall use the exact expressions for i/k, vlt, i^k^, 
and i/f^ in our investigation, we give also, for the purpose 
of comparison, the approximate expressions valid to first 
order of the dimensionless angular momentum j = J/M^, 
where J and M denotes the angular momentum and gravita- 
tional mass, respectively. The approximate expressions have 
the advantage to constrain the mass and the angular mo- 
mentum independently of the EOS. To first orde r in j the 
corresponding expression are (Miller et al., 1998b): 



1 

2^ 



M 

^circ 



M 



2_M^ 
2^rL~/ 



ms 
' circ 



ms 



:6M ( 1 - (f ) 2 j 



3 

6"2 1 

'2^M 

TT M' 



l + llx6-2j , 



(16) 

(17) 
(18) 

(19) 
(20) 



3 EQUATIONS OF STATE 
3.1 Neutron Star Matter 

The EOS of neutron star matter is the basic iiiput quantity 
to the structure equations discussed in Sect, bl Its knowl- 
edge over a wide range of densities is necessary. Whereas 
the density at the star's surface corresponds to the density 
of iron, the density in the centre of a very massive star can 
reach 15 times the density of normal nuclear matter. Since 
neutron star matter in chemical equilibrium (i.e. general- 
ized /3-equilibrium) is highly isospin-asymmetric and may 
carry net-strangeness its properties cannot be explored in 
laboratory experiments. Therefore, one is left with models 
for the EOS which depend on theoretically motivated as- 
sumptions and/or speculations about the behaviour of super 
dense matter. One main source of uncertainty is the com- 
petition between non-relativistic versus relativistic models. 
Though both treatments give satisfactory results for nor- 
mal nuclear densities, they provide quite different results if 
one extrapolates to higher densities (see, e.g, Ifuber et al. 
1998). Moreover, one encounters strong differences in the 
high density regime depending on the underlying dynamics, 
the many-body approximation and the assumptions about 
the composition. The simplest approach describes neutron 
star matter by pure neutron matter. Since pure neutron mat- 
ter is certainly not the true ground state of neutron star 
matter, it will quickly decay by means of the weak interac- 
tion into chemically equilibrated neutron star matter, whose 
fundamental constituents - besides neutrons - are protons. 



© 0000 RAS, MNRAS 000, 000-000 



4 Ch. Schaab and M. K. Weigel 



hyperons and possibly more massive baryons. Even a tran- 
sition to quark matter (so-called hybrid stars) and the oc- 
currence of pion- or kaon-condensates is possible. 

The cross section of a ne utron star can be split roughly 



into th ree dis tinct regimes (Borner, 1973; Weber & Glen 



dennin;;, 1993). The first one is the star's outer crust, which 
consists of a lattice of atomic nuclei and a Fermi liquid of 
relativistic, degenerate electrons. The inner crust extends 
from neutron drip density, p — 4.3 x 10^^ g c m~^, to a tran- 



sitio n de nsity of about pti — 1.7x 10 g cm' ( Pet hick et al, _ 
199£ |). l^ eyond this transition density pu one enters the star's 
third regime, that is, its core where all atomic nuclei have 
dissolved into their constituents, protons and neutrons. Fur- 
thermore, as outlined just above, due to the high Fermi pres- 
sure the core will contain hyperons, eventually more massive 
baryon resonances, and possibly a gas of free up, down and 
strange quarks. 

The EOS of the outer and inner crust has been studied 
in several investigations and is rather well known. We shall 
adopt for these regions the models derived by Ifaensel fc Pi- 



chon ( 1 )94 ) and by Negele fc Vautherin ( 1971: ), respectively. 
The models for the EO S of the star's core are discussed in 
detail in Schaab et al. (1996). 



An overview of the collection of EOSs used in this pa- 
per for the core region is given in Tab, h i . We have chosen 



a repre jentative collection of different models in order to 
check which EOSs are compatible with the QPOs and their 



theoret cal interpretation. 



An important characteristic of relativistic models is the 
appearance of a new saturation mechanism. In relativistic 
theories the repulsive force is caused by the exchange of 
vector mesons coupled to the baryon densities, whereas the 
attraction is coupled to the scalar densities by means of the 
scalar mesons. Hence, the repulsion increases with increasing 
density with respect to the attraction. Since non-relativistic 
treatments do not distinguish both kinds of densities, one 
has to introduce adhoc forces depending explicitly on density 
in order to reproduce properties of saturated nuclear matter. 

The potential parameters and coupling constants in the 
pure nucleonic sector are adjusted to nucleon-nucleon scat- 
tering data and properties of the deuteron in the case of 
the non-relativistic microscopic models and the relativistic 
Briickner-Hartree-Fock models. In this sense, these models 
are called parameter free. For the Hartree- and the Hartree- 
Fock approximation, as well as for the Thomas-Fermi model 
such an adjustment is not possible, since the free force pa- 
rameters would not yield saturation of nuclear matter. In 
such more phenomenological approximations the coupling 
constants are adjusted t o properties of saturat e d nuclear 



matter or atomic nuclei (Weber fc Weigel, 1989 



fc Mishnstin, 1996; Myers fc Swiatecki, 1996). In a more re - 
cent and more elaborate investigation ( Huber et al., 1998 ) 



Schaffnei 



the adjustment was performed with respect to properties 
of neutron star matter at 2-3 times nuclear density known 
from relativistic Briickner-Hartree-Fock calculations. Rela- 
tivistic Briickner-Hartree-Fock calculations cannot be per- 
formed over the whole density range at present, since an 
inclusion of hyperons leads to rather complex equation sys- 
tems. Even for pure nucleonic matter, the self consistent 
method is numerically stable up to 2-3 times nuclear den- 
sity only dHuber et al., 1998| ). 

Since the coupling constants of hyperons are not obtain- 



able from properties of normal nuclei and data of hypernu- 
clei are relatively scarce, one has a greater freedom in the 
selection of parameter sets in the hyperonic sector. From a 
theoretical standpoint one may first utilise the SU(6) sym- 
metry of the quark model for the vector coupling constants 
and adjust the cr-coupling from the binding energy of the 
A hyperon in nuclear matter. It turns out, however, that 
this procedure gives different cr-couplings depending on the 
chosen many-body approximation. Therefore, it seems rea- 
sonable to vary both couplings with the constraint of the 
A binding energy in nuclear matter in such phenomenolog- 
ical many-body theories. A further constraint in this pro- 
cedure is the compatibility with the allowed neutron star 
masses (for more details, see Huber et al. 1998 and Huber 
1998 ). The coupling constants for the strange mesons can 
be obtained to a certain extent from the d ata of double 
A-hyp ernuclei (Schaffner fc Mishustin, 1996). Batty et al. 



( 1994 ) claim some doubts on the existence of E hyperons in 
neutron star matter. However, the disappe arance of S hy- 



peron s would only slightly change the EOS ( Balberg & Gal 



1997). A further open question is still the existence of A iso- 
bars in the interior of neutron stars. In relativistic Hartree 
treatments, the isospin unfavoured A~ isobar does not ap- 
pear because of the rather large p-coupling, which is neces- 
sary to repr oduce the empirical symmetry coefficient ( Hubei 



et al., 1998). A isobars are therefore o ften neglected a prl 



ori in rel ativistic Hartree treatments ( schaffner fc Mishus- 



tin, 1996). However, in relativistic Hartree- Fock approxima- 



tions, the p-coupling becomes smaller due to the exchange 
contributions. For that reason, the charge favoured A~ may 
now enter the composition depending on the behaviour of 
the effecti ve A-mass in neutr on star matter, which is quite 
uncertain (Huber et al., 1995). 



The possibility of a transition of confined hadronic mat- 
ter into quark matter at high densities is included in the 
EOS labelled G^itS dGlendenning, 1995| ) (so-called hybrid 
stars). The transition was determined for a bag constant 
of B^^"^ — 180 MeV which places the energy per baryon of 
strange quark matter at 1100 MeV, well above the energy 
per nucleon in saturated nuclear matter as well as in the 
most stable nucleus, ^^Fe {E/A ^ 930 MeV). This model 
predicts a transition to a mixed phase consisting of quark 
matter and hadronic matter at a density as low as 1.6 tiq. 
The pure quark-matter phase is reached at a density of about 
10 no, which is close to the central density of the heaviest 
and thus most compact star constructed from such an EOS. 
A phase transition to pion condensation is accounted for in 
GJoo- This EOS predicts a phase transition at n ~ 1.3 tiq. 
The possibility of absolutely stable strange quark matter 
will be considered in the following section. 

The stiffness of the EOS depends strongly on the in- 
ternal degrees of freedom. Generally, the more degrees of 
freedom are taken into account the softer the EOS becomes 
(s. Fig. ^. A softer EOS, in turn, leads to lower maximum 
neutron star masses and, for fixed mass, to higher central 
densities. 



3.2 Strange Stars 

The hypothesis that strange quark matter may be the ab- 
solute ground state of the strong interac tion ( not ^^Fe) has 
been raised independently by Bodmer (1971) and Witten 
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Table 1. Dynamics and approximation schemes for EOSs derived for the cores of neutron stars 



EOS 



Composition 



Interaction 



Many body approacli 



UV14 + UVII 
UV14+TNI 
TP 
HVpn 
HV 
RHl 
RHFp„ 
RHPl 
RHFS 
G300 

r.i{240 
^M78 
f-iK240 

r.K240 

^B180 
TMl-ml 
TMl-m2 

NLSHl 

NLSH2 
SMb145 
SMb160 



P, n, 
A, 

A, 

A, E^ 
p, n. 



p, n, e , ^ 
p, n, c~ , ^~ 
p, n, e~ , ^~ 
p, n, e~ , ^~ 
A, E±.0, =0--, e-, ^ 

p, n, , ^~ 
A,' E±.o', h".", e-, f 



p, n 


A, E±'0, 


gO,- 


e 




p, n 


A, E±.o 


gO,- 


e~ 




, A, 




~ 1 o~ 


M~ 


■ q 


p, n 


A, e'±.0, 


gO,- 






p, n 


A, E±.0, 


gO,- 






p, n 


A, E±'0 


=0,- 


c 




p, n 


A, E±.0 


=0,- 


c 





7r-condcnsation 



, quarks 



Urbana V14 and Urbana VII 
Urbana V14 and TNI 
TF96 
exchange of a, uj , p 
exchange of (j, uj , p 
exchange of cr, 7r, p 
exchange of cr, 7r, p 
exchange of tj, Ci), 7r, p 
exchange of tj, c^;, 7r, p 
exchange of (j, lj, p 
exchange of (j, lj, p 
exchange of (j, oj, p 
exchange of (j, uj , p 
exchange of (j, lj, p 
exchange of (j, oj, p 
cxclrangc of cr, oj, p, 
exchange of (7, p, (j*, 

exchange of (j, oj, p, 
exchange of ct, CiJ, p, (j , 
Bag model 
Bag model 



NRV 
NRV 
Thomas-Fermi mo 
RH 
RH 
RBHF+RH 
RBHF+RH 
RBHF+RHF 
RBHF+RHF 
RH 
RH 
RH 
RH 
RH 
RH 
RH 
RH 
RH 
RH 
Fermi gas 
Fermi gas 



Wiringa, Fiks & Fabr' 

Wiringa et al. t|l9SHl 



(LStki 



1 QttQ 



1 QQr^ 



1997 



Strobel et al. ( ipS*] 
Weber &: Weigel 
Weber &: Wcige 
Huber ct al 
Huber ct al 
Huber ct al 
Huber et al. r \aai 
Glendenning 
Glendcnning 
Glendenning (1^^ 
Glendcnning 
Glendenning 
Glendenning 
Schaffner & Mishustin (liaSd) 
Schaffner & Mishustin ( E55I ) 
Schaffner & Mishustin ( E551 ) 
Schaffner & Mishu qtin ( 199e ) 



Farhi & Jaffe ( M»4 ) 
Farhi & Jaffe (1984) 



Abbreviations: NRV=non-relativistic variational method, RH=relativistic Hartree approximation, 
Hartree-Fock approximation, RBHF=rclativistic Briickncr-Hartree-Fock approximation. 



RHF=relativistic 




e [g cm " 

Figure 1. Pressure-density relation for several EOSs 



(1984). If the hypothesis is true, then a separate class of 
compact stars could exist, which range from dense strange 



star s to strange dwarfs to strange planets ( Weber et al. 
1995; G lendenning et al., 1995b; Glendenning et al., 1995a) 



In principle both strange and neutron stars could exist. How- 
ever if strange stars exist, the galaxy is likely to be contami- 
nated by strange quark nuggets which would convert via im- 



pact "conventional" neutron stars into st range stars (Glen- 



dennin;;, 199C; Madsen & Olesen, 1991; Caldwell & Fried 



man, 1! '91). This would imply that the objects known to as- 



tronomers as pulsars are probably rotating strange matter 
stars and not neutron matter stars, as it is usually assumed. 
Presently there is no sound scientific basis on which one can 
either confirm or reject the hypothesis, so that it remains a 
seriou s possibility of fundamental significance for compact 
stars ( Weber et al., 1997b ; Weber et al., 1997a ; Madsen 



1990 plus va rious other physical phenomena (Madsen & 
Hacnsel, 1991). Below we will explore the implications of 



the existence of strange stars with respect to the QPO phe- 
nomena. This enables one to test the possible existence of 
strange stars and thus draw definitive conclusions about the 
true ground state of strongly interacti ng m atter. 

As pointed out by Alcock et al. (1986), a strange star 
can carry a solid nuclear crust whose density is strictly lim- 
ited by neutron drip. This possibility is caused by the dis- 
placement of electrons at the surface of the strange matter 
core, which generates an electric dipole layer there. It can 
be sufficiently strong to produce a gap between ordinary 
atomic matter (crust) and the quark-matter surface, which 
prevents a conversion of ordinary atomic matter into the 
assumed lower-lying ground state of strange matter. Obvi- 
ously, free neutrons, being electrically charge neutral, can- 
not exist in the crust, because they do not feel the Coulomb 
barrier and thus would gravitate toward the strange-quark 
matter core, where they are converted by hypothesis into 
strange matter. Consequently, the density at the base of the 
crust (inner crust density) must always be smaller than neu- 
tron drip density. The situation is illustrated in Fig. |^ where 
the EOS of a strange star with crust is shown. 

The strange-star models presented in the following are 
construc ted fo r an EOS of strange matter derived by Farhi 
& Jaffe (1984). We shall study models with a fixed strange 
quark mass, rUs = 100 MeV, and the two bag constants, 
B^/-* = 145 MeV and B^^^ = 158 MeV. 



4 KILOHERTZ QUASI PERIODIC 
OSCILLATIONS 

In Table ^ the observations of kilohertz QPOs from 16 
sources are summarised. The sources can be classified into 
two classes, the Atoll- and the Z-sources, depending on the 
form of their colour-colour diagram. Almost all sources show 
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Table 2. Observational data of kilohcrtz QPOs 



oourcc 


lype 




ru 1 


A ru 1 


t^Burst L-H-ZJ 


JrvcicrcncGS 


4U 0614+091 


Atoll 


750-1145 ± 10 


480-800 


323 ±4 


328t 


1 


4U 1608-521 


Atoll 


940-1125 


650-890 


233 ± 12-293 ± 7 




2,3,4,5,6 


4U 1636-5361 


Atoll 


1147-1228 


830-940 


257 ± 20-276 ± 10 


581 


6,7,8 


4U 1728-34 


Atoll 


988-1058 ± 12 


638-716 


355 ±5 


363 


7,9 


KS 1731-260 


Atoll 


1176-1197 ± 10 


900 


260 ± 10 


523.92 ± 0.05 


10,11 


4U 1735-44 


Atoll 


982-1161 ± 1 


632-729 


296 ± 12-341 ± 7 




12,13 


4U 1820-30t 


Atoll 


640-1060 ± 20 


400-800 


275 ±8 




14 


Aql X-1 


Atoll 


( oou 








1 ^ 


GX 3-1-1 


Atoll 










16 


4U 1705-44 


Atoll 


1074 ± 10 


776-867 


298 ± 11 




17 


Sco X-1 


Z 


872-1115 


565-890 


310-230 




18 


GX 5-1 


Z 


567-895 


325-448 


298 ± 11 




19 


GX 17-1-2 


z 


645-1087 ± 12 


480-781 


294 ±8 




20 


Cyg X-2 


z 


731-1007 ± 12 


490-780 


346 ± 29 




21 


GX 349-1-2 


z 


978 ±9 


712 


266 ± 13 




22 


GX 340-f 


z 


567-820 ± 19 


247-625 


325 ± 10 




23 



& Bradt ( 1997 ), _L2LWijnands et al. (|l998d), 13: Ford ct 
: 6: Str ohmaver (1998), 17: Fo rd, va n der K 



al. (199i 



f: Q PO is statistical not significant, |: i^qpoi is ind epende nt of luminosit y. Refe rences 1: Ford et a l. (|1997| ), 2: Mendez et 
al. (|l998a| ), 3: Mendez et al. t|l998h| ), 4: Yu et al. ( |l997a|). 5: Yu ct al. ( 1997b| ), 6: Kaaret et al. (|l997^ 7: van der Klis 
(|l99^). 8: _Stroh mayer et al. (|l998|). 9: Stroh mayer et al. (|l99 d). 10: Wijnands & van d er Klis (1997), 11: Smith, Morga n 

|998h), 14: Zhang et al. ( 1998b| ) 15: Zha,ng et al. ( 1998a| ), 
Kaaret ( |l998a|) 18: van der Klis & Wijnands ( |l997|. 19: Wiinands et al. 
( |l998ht),_ 20: Wijnands et al. ( |1997| ), 21: Wijnands et al. ( |l998a ), 22: Zhang, Strohmayer & Swank ( P98c ), 23: Jonker et 
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Figure 2. Pressure-density relation for strange stars with crust. 
The bag constant is B^/'' = 145 MeV or 158 MeV, respectively, 
the mass of the strange quark is rris = 100 MeV. 



simultaneously two kilohertz QPOs. Although the frequen- 
cies of both QPOs vary by several hundred Hertz with the 
accretion rate, the variation of the frequency separation is 
only small. 

Therefore the observations strongly favour some kind 
of beat-frequency model. In such models the upper kilohertz 
QPO corresponds to the Kepler rotation at the inner edge of 
the accretion disk. The lower QPO corresponds to the beat 
frequency between the upper QPO and the spin of the neu- 
tron star. This interpretation is supported by the fact that 
the frequency separation is equal to the QPO frequency (or 
half of it) in type I X-ray bursts observed in some of the Atoll 



sources (see Tab. The question remains how ever w here 
the upper QPO is generated. Strohmayer et al. ( 1996 ) sug- 
gested that this radius is id entical to the magnetospheric 
radius, whereas Miller et al. (1998c) propose that it is iden- 
tical to the sonic radius (see van der Klis 1997 for a critical 
discussion of these models). In both cases, the orbital radius 
has to be l arge r than the radius of the innermost stable orbit 
(see Sect. 2.2), or equivalently, the frequency vqpoi of the 
upper QPO has to be smaller than the Kepler frequency at 
the innermost stable orbit: 



I'QPOi < ~ 2200^ Hz, 



(21) 



where the approximate expression (^^ was used. The exact 
expression depends on the EOS and on the spin period of the 
neutron star. This inequality sets an upper limit to the mass 
of the star. Obviously, the orbital radius has to be larger 
than the star's radius, which is larger than the radius of the 
innermost stable orbit for low star masses. This constraint 
sets a lower limit to the star mass which again depends on 
the spin period and the EOS. 

Figure ^ shows the maximally allowed Kepler frequency 
as function of the gravitational mass for two EOSs, RHFpn 
and RHF8. If the radius of the innermost stable orbit is 
smaller than the star's radius the maximally allowed Kepler 
frequency is given by the Kepler frequency at the star's equa- 
tor. This kind of solution is represented by the rising branch 
of UKiM), since increases with increasing gravitational 
mass. If the innermost stable orbit is located outside the 
star the maximally allowed Kepler frequency decreases again 
with increasing mass. The two EOSs differ for M > 1.0 M© 
in the composition. In RHF8 hyperons are included, whereas 
RHFpn assumes pure nucleonic matter. RHF8 is therefore 
softer than RHFpn at high densities. This means also, that 



© 0000 RAS, MNRAS 000, 000-000 



Quasi Periodic Oscillations and Constraints on the Equation of State 7 



1600 



1400 - 



1200 



1000 



800 - 




600 



0.5 



1.0 



1.5 
M[M3, 



2.0 



2.5 



Figure 3. Maximally allowed Kepler frequency ^ function of 
the neutron star's mass M for the two EOSs RHFpn and RHF8 
and two different spin frequencies = 275 and 355 Hz. 



RHF8 can only support smaller masses than RHFpn. The 
maximally allowed Kepler frequency depends only slightly 
on the neutron star's spin frequency, as long as the spin fre- 
quency is not larger than about ten percent of the star's 
Kepler frequency. If one compares the !/K-curves for exam- 
ple with the frequency i'qpoi = 1060 of the highest kilo- 
hertz QPO observed in 4U 1820-30 (;^s = 275 Hz), one ob- 
tains that the mass of the neutron star has to be within the 
range between M = 1.0 M© and 2.25 M© (for RHFpn) and 
Mniax = 1.65 Mq (for RHF8), respectively. 

If one adopts the interpretation of the sonic-point model 
one expects that the inner boundary of the accretion disc is 
close to the innermost stable orbit. Indeed, the observations 
of the sources 4U 1820-30, 4U 1608-52, and 4U 1636-536 
seem to confirm this interpretation, since the Q PO frequency 



rem ains constant for high m ass accretion rates (Kaaret et al. 
1997|; [Z|hang et al., 1998h| ). Therefore, only the right inter- 



section point of the i/K(Af)-curve with the hue uk ~ t'qpoi 
is allowed by the observations|||. The mass of the neutron star 
in the binary system 4U 1820-30 can now be determined. For 
RHFpn the mass is equal to M = 2.25Mo whereas RHF8 
would be excluded. 

The figures ^-|^ show the respective ranges of masses 
for which the Kepler frequency at the equator or at the 
innermost stable orbit is larger than the highest observed 

8PO frequency i'qpoi for a-U EOSs considered here (see Tab. 
. The vertical line refers to the approximate mass, which 
is obtained by setting j = in Eq. ([l£|). This expression 
underestimated the upper limit of the mass. The dashed 
bars refer to the EOSs for which models with i/^" = i^qpoi 
do not exist, i.e. whose maximally stable mass is to small. If 
the interpretation of the highest observed QPO frequency in 
the three sources 4U 1820-30, 4U 1608-52, and 4U 1636-536 
is confirmed, the mass of the neutron star is constrained to 



t The left intersection point is unlikely, since it is expected that 
the lower and the upper frequency QPO cannot be observed at 
the same time if the accretion disk would extend to the neutron 



star surface (Zhang et al., 1997). This would be contrary to the 
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Figure 4. Range of masses for which the Kepler frequency at 
the equator or at the innermost stable orbit is larger than the 
highest observed QPO frequency i^qpoi ~ 1228 Hz in the source 
4U 1636-536 {us = 290 Hz). The dashed bars refer to the EOSs for 

do not exist. The approximate 



value Ma 



QPOi 

I.SOMq is also shown. 



the respective right end of the solid bars. The EOSs, for 
which models with = '^qpoi <io not exist (dashed bars) 
are excluded. 

Especially, the observation of 4U 1820-30 agrees only 
with the two nucleonic EOSs UV14-I-UVII and RHFpn. If 
the interpretation of the highest observed QPO frequency 
in this source is confirmed by further observations, the ex- 
istence of hyperons, meson condensates, quark condensates, 
and strange stars can be rejected. The neutron star mat- 
ter consists thus only of nucleons and leptons. At least, the 
hyperon fraction (i.e. strangeness) is much smaller than pre- 
dicted by actual relativistic calculations. 

Figures ^ and ^ show the allowed ranges of mass of 
the neutron stars in the respective binaries for the EOSs 
RHFpn and RHF8, respectively. For RHF8, the masses are 
limited by the maximally stable star mass Mmax ~ 1.65 Mq, 
which depends on the spin period. All sources are consistent 
with a canonical star mass of M = 1.4Mq. Only if one 
assumes that the highest QPO frequency is equal to the 
Kepler frequency at the innermost stable orbit, the masses 
lie between M = 2.0 Mq and 2.5 Mq (for RHFpn). 



CONSTRAINTS ON EQUATION OF STATE 
BY LENSE-THIRRING PRECESSION 



simultanous observation of both QPOs in several sources. 



Stella & Vietri ( |199^ ) suggested that a third oscillation with 
a frequency vqpos around ~ 10 Hz is also produced at the 
inner border of the accretion disk by Lense-Thirrine pre- 
cession. The observed frequencies i'qpos (s. Tab. ra) can 
again be compared with theoretical neutron star models. 
The neutron star models are constructed for a given spin 
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Table 3. Observational data of QPOs in the 10 Hz range. 



Source Vs [Hz] uyl [Hz] i^QP03 [HZ] 



References 



4U 0641+91 323 



900 



22 



Stella & Vietri (1998) 



4U 1608-52 233 



800 



20 



Yu et al. (1997E) 



4U 1728-34 355 



355 
551 
699 
1122 



9.0 
14.1 
26.5 
41.5 



Ford & van der Klis (1998) 



KS 1731-260 



262 



1197 



27 



Wijnands & van der Klis (1997) 



4U 1735-44 



326 



1149 



29 



Wijnands et al. (1998c 



Sco X-1 



247 



1050 
1101 
1135 



5 
9 

13 



van der Klis (1996) 



GX 17+2 



294 



644 
832 
1042 



5 
9 

13 



Wijnands et al. (1997) 



Cyg X-2 



346 



731 
856 
1007 



5 
9 

13 



Wijnands et al. (1998a) 



GX 340+0 



325 



569 
730 
823 



5 
9 

13 



Jonker et al. ( 199^ ) 



GX 5-1 



298 



506 
685 
889 



5 
9 

13 



Wijnands et al. ( |l998b| ) 



UV„+UVII 
UV„+TNI 
TF 
HV„ 

pn 

HV 
RH1 

RHFp„ 
RHF1 
RHF8 

^300 

TM1-m1 
TM1-m2 
NLSH1 
NLSH2 

h 

h 



K240 
K240 

B 

K240 



I 

I-- 



0.0 



0.5 



H 



I 1 



1 

h 

I 

h 

I 



1.0 



1.5 



2.0 



M[M3, 



2.5 



UV,.+UVII I- 

UV,.+TNI V 

TF 

HVp„ 

HV 

RH1 

RHFp„ 

RHF1 

RHF8 

^300 
_ K240 

TM1-m1 
TM1-m2 
NLSH1 
NLSH2 

I- 

h 



K240 
K240 



0.0 



0.5 



1.0 



1.5 



2.0 



M[M3, 



2.5 



Figure 5. Same as in Fig. i] but for the source 4U 1728-34 



max 

QPOl 



1100 Hz, Vs = 355 H 



A/a 



: 2.01Mq). 



Figure 6. Same as in Fig. ^ but for the source 4U 1820-30 

'..max ----- 
'^QPOl 



frequency Us and a fixed, unfortunately unknown, star mass 
M. The obtained monotone relations VY^{r,6 = 7r/2) and 
vur{r,6 — 7r/2) are transformed into tiie relation vurii^K)- 



Since tlie Lense-Thirring precession frequency ult is, in first 
approximation, proportional to i^s tiie relation uut/i^b as 
function of vk is shown in Fig. S for all EOSs. The observa- 
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Figure 7. Same as in Fig. |^ but for different sources and the 
fixed EOS RHFnn. 
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Figure 8. Same as in Fig. I but for RHF8. 



tions are shown as circles. The curves depend only weakly 
on the neutron star mass. One can therefore draw conclu- 
sions from the comparison of the theoretical curves with the 
observations, although the star masses are unknown. Since 
both virsjvs and vy^ do not depend on in first approxima- 
tion, the calculation of the theoretical curves for one specific 
spin frequency, — 363 Hz (this corresponds to the spin 
frequency of the neutron star in 4U1728-34), is sufficient. 

As can be seen in Fig. ^, the observations of KS 1731-260 
and 4U 1735-44 agree with the theoretical curves of group 
1, which contains the stiffest EOSs of our sample. This is in 



agreement with the results of Stella & Vietri ( |l998 ). How- 
ever, the other observations, including the observations of 
the Z-sources which are not shown, lie above all theoretical 
curves. The observed frequency i^Qpos is thus higher than 
the Lense-Thirring precession frequency even for the stiffest 
EOSs. If one assume that the observed frequencies corre- 
spond to the first overtone of the Lense-Thirring precession, 
the Atoll-source observations are in the range of the theo- 
retical curves. 

In the case of the Z-sources, the detected frequencies 
J'QPOS are not only too large in most sources, but addition- 
ally the slope of the relation i'Qpo3(i'QPOi) of the Z-source 
observations is much higher than the slope of ^'lt(^'k) for the 
theoretically determined models. Even the assumption, that 
the first overton e of the Lense-Thirrin g precession frequency 
i^LT is detected (Stella & Vietri, 1998) cannot explain these 
discrepancies. 



6 CONCLUSIONS 

In this paper, we derived models of rapidly rotating neutron 
stars and strange stars by solving the general relativistic 
structure equations for a broad collection of modern EOSs. 
We compared the space time geometry of these models with 
recently discovered QPOs in the X-ray brightness of LMXBs. 

If one follows the general beat-frequency interpretation 
of the kilohertz-QPOs, i.e. that the higher frequency QPO 
originates at a stable circular orbit, one can constrain the 
mass of the neutron star to a range which depends on the 
EOS. This mass range is for all sources and for all EOSs 
consistent with a canonical mass, M — 1.4 Mq. As it was 



stated by Miller et al. ( [1998a| ) and Thampan et al. ( |l99£| ) 
the exact lower and upper limits of the neutron star's mass 
can only be determined by using fully relativistic models of 
rapidly rotating neutron stars. The exact limits differ from 
the approximations with j = by roughl y 10 % ■ 

As it wa s shown by Kaaret et al. (1997) and Zhang 
et al. (1998b), the observation of a maximum frequency 
i/QPOi of the high frequency QPO in the sources 4U 1820- 
30, 4U1608-52, and 4U 1636-536 favour the interpretation 
that this QPO originates at the innermost stable orbit. 
If this interpretation is correct, the mass of the neutron 
star can be exactly (within the observational errors) deter- 
mined for a given EOS. The approximately obtained mass, 
M = 2.07 Mq, of the source 4U 1820-30 is larger than the 
maximum mass of most of the considered EOSs. This con- 
clusion is even strengthened if the observed maximum fre- 
quency i^Qpoi is compared with the exact neutron star mod- 
els. The only allowed EOSs of our broad collection are then 
the UVi4-f-UVII and RHFpn, which both describe neutron 
star matter as consisting of nucleons and leptons only. The 
derived masses of the three sources lie in the narrow range 
between M = 1.92 and 2.25 Mq. 

This result is also of some importance for the nature 
of the object left in the supernova 1987A. During the first 
t en seconds afte r the supernova, neutrinos were detected 
([ Chevafier, l"997| ). This means that a protoneutron star was 



formed in the supernova. The fact that up to now no pul- 
sar emissi on co uld be detected was interpreted by Bethe 
& Brown (1995) that the protoneutron star collapsed to a 
black hole when the star became transparent to neutrinos 
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Figure 9. Ratio vuy/vs of the Lense-Thirring precession frequency !/lt ant the spin frequency = 363 Hz as function of the Kepler 
frequency u^. The used mass is M = IAMq. The first group of curves contains the EOSs: NLSHl, NLSH2, TMl-m2, TMl-ml, RHFpn, 
G30O1 HVpn, and HV (from the left to the right; the bold typed label corresponds to the bold curve). Group 2: Gp^'"', RHF8, RHl, 
G30O' GmtIu: and GK2|°. Group 3: SMbi45, TF, UV14+UVII, UV14+TNL and G|2|0. Group 4: RHFl and SMsieo- The full (open) 
circles correspond to (half of) the observations of the Atoll-sources in Tab. til 



after roughly 10 s. The estimated value of the baryonic mass 
Mb ~ 1.63 - 1.76 Aff^ (iBethe fc Brown, 1995|; iThiclcmann 



probing the interior of neutron stars. Compare d to the other 
tools like o bservations of the maximu m mass ( yan Kerkwiik 



et al. , 1 39€ ) gives thus an upper limit to the maximum grav- 



el al., 199i:), the limiting spin period (Friedman et al., 198(: 



is in cleir contradiction to the derived mass of, e.g. 4U 1820- 
30. 

It is generally believed that neutron stars are born 
with a mass about 1.4 — 1.5 M©. Neutron stars in LMXBs 
would therefore accrete 0.4 — 0.8 Mq during their lifetime. 
It seems reasonable to assume that some of the neutron 
stars in LMXBs accrete enough matter to reach the max- 
imally stable mass (Mmax = 2.20 for UV14+UVII, 
Mmax = 2.44 M© for RHFpn). The neutron star would then 
collapse to a black hole. 

The interpretation of the QPO with frequencies t'qpos 
about 10 Hz as Lense-Thirring precession frequency of the 
accretion disk seems not to be consistent with the theoretical 
star models, unless one assumes that the first overtone of the 
precession is observed. In the case of Z-sources however, the 
necessary ratio of the frame dragging frequency and spin 
frequency, and thus the moment of inertia, of the models if 
too small compared to the observed frequency i^qpos or half 
of it. 

Our general conclusion is that the observations of kilo- 
hertz QPOs in LMXBs provide us another powerful tool for 



itation al mass of a neutron star: Mmax 5, 1-6 Mp). This limit Weber fc Glendenning, 19921) , an d cooling sirn ulations ( Fsu 



ruta, 196C; Schaab et al., 1996| ; Page, 1998), the derived 



constraints, especially in the sonic-point-interpretation, are 
rather strong. The lower limit Mmax ^ 2.15 Mq is only con- 
sistent with two EOSs, UV14-I-UVII and RHFpn, which are 
relatively stiff at high densities. Their stiff behaviour at high 
densities seems to be only possible if the neutron star matter 
consists of neutrons, protons, and leptons only. At the most, 
a small admixture of hyperons may be allowed. However, one 
has to admit that such a composition somehow contradicts 
our conception of neutron star matter, since fieldtheoreti- 
cal models lead more or less inevitably to a more complex 
composition at high density. If the given interpretation is 
correct, one has therefore to reinvestigate the problem of 
super dense neutron star matter by dropping, for instance, 
the standard assumptions about the coupling of the hyper- 
ons in such regimes. 

This specific conclusion depends however on the inter- 
pretation of the majcimum frequency of the kilohertz QPO 
within the sonic-point-model. As a kind of beat-frequency 
model, this model predicts a constant frequency separa- 
tion = !^Qpoi — i^QP02- Recently two further exam- 
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pies, where this constancy is not ob served, has been dis- 



covered: 4U 1608-5 2 and 4U 1735-44 (Mendez et al., 1998c 



Ford et al., 1998b). Moreover, the frequency separation in 
the Atoll source 4U 1636-536 seems not to be consistent 
with the half of the frequency of the QPO in type I bursts 
(Mendez & van Paradijs, 1998). Though the beat-frequency 



models are the most hopeful candidates in explaining the 
QPO-phenomenology it has to be awaited how the variation 
of the frequency separation and deviation from the QPO 
frequency in bursts can be incorporated to these models. 
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